You can find all related materials, links and references there.

• Missing data is ubiquitous (Allison 2001).
• Ad hoc solutions may yield invalid inferences (Van Buuren 2018).
• Rubin (1987) proposed the framework of MI.

simulation <- function() {
  mvtnorm(X,Z1,Z2) %>% 
  mutate(Y~X+Z1+Z2) %>% 
  for (max_iterations in 1:100) {
    ampute() %>% 
    impute() %T% 
    convergence_diagnostics %>% 
    lm(Y~X+Z1+Z2) %>% 
    pool %>% 
    simulation_diagnostics %>% 
    c(., convergence_diagnostics)
  }
}

replicate(simulation, n = 1000)

[]